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ABSTRACT 

I present models of disk evolution around young binary stars. I show that the primary factor in determining 
circumbinary disk lifetimes is the rate of disk photoevaporation. I also find that photoevaporative clearing 
leaves a signature on the distribution of circumbinary disk lifetimes, with a sharp increase in disk lifetimes for 
binary separations a <0.3— 1AU. Observations of young binary stars can therefore be used to test models of 
disk evolution, and I show that current data set a strong upper limit to the rate of on-going photoevaporation 
(< lO~ 9 M0yr -1 ). Finally I discuss the implications of these results for planet formation, and suggest that 
circumbinary planets around close {a < 1 AU) binaries should be relatively common. 

Subject headings: planetary systems - binaries: close - protoplanetary disks - planets and satellites: formation 
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1. INTRODUCTION 

It has long been recognised that mos t Sun-like stars form 
in bi nary or multiple systems (e.g., iDuquennoy & Mayor 
1199 lb . Circumbinary disks are a natural con sequence of 
binary star formation (e.g., iMonin et alj 120071) . and these 
young disks are potentially the sites of planet formation. 
Recently the Kepler mission has disco v ered a number of 
circumbinary plan ets dDovle et alj 1201 U IWelsh et"aTJ 120121 ; 
lOrosz et al]|2012allbl) . the existence of which poses interest- 
ing challenges for planet formation theories. Moreover, new 
high-resolution techniques have now begun to build a useful 
censu s of disks around, young, ~AU-sepa ration binary stars 
(e.g.. iKraus et alj|2012t Harris et al] 120121) . These observa- 
tions suggest that planets form readily in circumbinary disks, 
but little is known about the physical conditions in these com- 
plex young systems. 

In this Letter I present a simple one-dimensional (1-D) 
model for the evolution of protoplanetary disks around close 
(< 10AU) binary stars. The formation and early evolu- 
tion of such disks is dominated by gravita tional instabili- 
ties and magnetically-driv en outflows (e.g., iDurisenl 1201 li 
iKonigl & Salmeronl 1201 ll) . while later evolution and final 
disk dispersal is driven by the competiti on between vis- 
cous accretion and pho toevaporation (e.g. JClarke et alj|2001t 
I Alexander et al.ll2006bl) . However, the tidal torque from a bi- 
nary inhibits disk accretion and modifies the evolution sub- 
stantially from the single-star case. I find that the rate of 
disk photoevaporation plays a critical role in determining how 
young circumbinary disks evolve, and that photoevaporative 
clearing leaves a characteristic signature on the distribution 
of disk lifetimes. I also show how observations of circumbi- 
nary disks can be used to inform our understanding of proto- 
planetary disk clearing, and discuss the consequences of these 
results for planet formation around binary stars. 

2. MODEL 

The evolution of an accretion disk ar ound a binary is de- 
scribed by (e.g., Lin & Papaloizou 1986) 
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Here £(/?, t) is the disk surface density, t is time, R is the cylin- 
drical radius (relative to the centre of mass), v is the kinematic 
viscosity, M tot is the total stellar mass (i.e., M tot = M\ +M2, 
where M\ & M2 are the masses of the individual stars), and 
S w (/?,f) is the mass-loss due to photoevaporation. The binary 
orbit is assumed to be circular and co-planar with the disk. 
For a binary of mass ratio q = M2/M1 and semi-major axis a, 
the ti dal torque from the b inary, A(R,a), is approximated by 1 
(e.g.. lArmitage et al.ll2002l) 



A(R,a) = 



q 2 GM tot ( a 
2R 



where 



Ap = max(H, \R-a\) 



(2) 
(3) 



and H is the disk scale-height. The back-reaction of this 
torque causes the binary orbit to shrink, but this depends sen- 
sitively on the gas dynamics inside the inner edge of the disk 
(which cannot be accurately modeled in 1-D). Moreover, the 
bulk of this evolution happens at early times, while here our 
primary interest is the late stages of the disk evolution, so for 
simplicity I fix the binary separation a to be constant through- 
out. Test calculations show that this approximation does not 
influence the results significantly. 

Accretion in protoplanetary disks is thoug ht to be drive n 
by magnetohydrodynamic turbulence (e.g., iBalbusi 1201 ll) . 
Here this is approximated by an alpha-prescription for 
the disk viscosity v{R) = a£lH 2 , where a = 0.01 is the 
IShakur a & Sun yaevl (1 19731) viscosity parameter and Vt{R) = 
y/ GM tot /R 3 . a here essentially represents a time and space 
average of the efficiency of angular momentum transport in 
the disk, and the choice a = 0.01 is consiste nt with observa- 
tions of protoplanetary disk accretion (e.g., lHartmann et al.l 
1998). I choose a power-law form for the disk scale-height 
H oc R 5 / 4 , which results in a linear viscosity law v oc R. The 
power-law is normalized by setting H/R = 0.05 at R = 1 AU. 

It is now well established that photoevaporation by high- 
energy phot ons dominates protoplanet a ry disk clearing at lat e 
times (e.g., iPascucci & Sterzikl l2009t iPascucci et al.l [20111) . 
This radiation heats the disk surface, and beyond some crit- 
ical radius the heated gas is unbound and flows as a wind. 

1 Strictly this expression for A applies only for R> a, but I assume E = 
interior to the binary orbit so the case R < a is not relevant here. 
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When the disk is optically thick to the high-energy radia- 
tion this mass-loss is concentrated around the critical radius 
R c ~ 0.2GMt ot /c 2 , where c s is the sound speed of the heated 
disk atmosphere. However, if there is an optically thin inner 
cavity the inner disk edge can be photoevaporated directly, 
and the mass-loss is instead concentrated close to the disk 
edge. Around single stars this occurs only at late times, when 
the wind is able to overcome disk accretion. Here, however, 
the torque from the binary can clear the inner disk, and change 
the qualitative behaviour of the photoevaporative wind. 

Unfortunately i t is not yet clear whet her X-ray (e.g., 
Owen et ail 120101) or UV irradiation (e.g., Alexander et al. 
2006b) drives the photoevaporative wind. However, the mass- 
loss profile is similar in both cases, so for simplicity I adopt 
a parametrized form for the photoeva porative mass-loss term 
(motivated by analyti cal models, e.g., Hollenbach et al.lll994t 
lAlexander et al.l2006al) . When the inner disk is optically thick 

and when the inner disk is optically thin (i.e., the disk inner 
edge lies at /?,„ > R c ), we instead have 

With this form the mass-loss per unit area peaks at R c (in 
the optically thick case), and M t hick and M tn in are the inte- 
grated mass-loss rates (normalised to R m = 2R C in the optically 
thin case). The second form for £(/?) is used when the sur- 
face density interior to R c falls below the critical value EtMck, 
and this parametrization successfully mimics the behavior of 
more sophisticated models. I define two wind models: one 
with a low wind rate and small critical radius, and a sec- 
ond with a higher wind rate and a larger critical radius. The 
first represents photoevaporation by ioni zing (EUV) photons 
dFont etal.ll2004t lAlexander e"t~ai1l2006al) . with R c = 1.4AU, 
Mthick = 1.6 x lO-^Mgyr- 1 , = 1.1 x lO^Mgyr^and 

StMck = 10~ 5 gc m~ 2 . The second repre sents photoevapora- 
tion by X-rays (lOwen etal.ll20Tol l201lh . with R c = 5.0AU, 
Mthick = Mthin = 1.0 x lO^Mgyr" 1 , and E thick = lO^gcm" 2 . 
These two wind models are henceforth referred to as the 
"weak" and "strong" photoevaporative winds, respectively. 

The form of the torque function A does not allow any gas to 
accrete on to the binary, and in the absenc e of mass-loss the 
model therefore describes a decretion disk (Pringle Tl991l) . In 
practice, however, we expect gas to accrete from the disk on 
to the binary via tidal streams. This accretion flow is variable 
and is modulated by the binary orbit, but the average accretion 
rate on to the binary can be as much as ~ 10% of the steady- 
state disk accretion rate (e.g., iMacFadyen & Milosavrjevia 
2008). I therefore allow gas to accrete from the inner disk 
edge on to the binary at a fraction e of the disk accretion rate 
(computed as 37rz/£ at R = 3R m ). These accretion streams 
may be optically thick to the radiation that drives photoevap- 
oration, but detailed modeling of how the streams "shield" the 
wind is beyond the scope of this initial investigation. Instead 
I define an average surface density S streams = 2M; n /Slb/? 2 n , 
where M m is the accretion rate from the inner disk edge and 
fib is the orbital frequency of the binary, and switch wind pro- 
files only when E streams < S thic k. 

Here we are primarily interested in the late-time evolution 
of the disk, so I adopt a simplified set of initial conditions. The 



initial disk mass is taken to be Md = 10 _1 5 M©, and the ini- 
tial surface density profile is assumed to be an exponen tially - 
truncated power-law (e.g.. lLvnden-Bell & Pringlell 19741) 

E( * )= ^ exp (-£) • (6) 

The inner disk edge is initially a step-function at 5a, but 
rapidly relaxes to a self-consistent profile. The scaling radius 
Ru determines the initial disk size and is set, arbitrarily, to be 
Ro = 15a (i.e., three times the inner edge radius). This implic- 
itly assumes that the disk angular momentum scales with that 
of the binary; this assumption is not well justified, but in the 
absence of a first-principles model for binary formation this is 
the most sensible way to proceed. 

This set of disk models is therefore described by four pa- 
rameters: the binary separation a and mass ratio q, the pho- 
toevaporation model (strong or weak), and the accretion effi- 
ciency e. I have run grids of models with values of log I0 (a) = 
-1.0, -0.99, -0.98. ..1.5 (i.e., a = 0.1-31.6AU). The standard 
model grid uses an equal mass binary (q = 1) and allows gas 
to accrete on to the binary with e = 0.1. Variant models were 
also run with q = 0.3, e = and e = 0.01, as well as a reference 
set of single-star models (which have A = and the same Ro 
as the standard models). All of the model grids were run with 
both the strong and weak photoevaporative winds. The mod- 
els were integrated forwards in time until the circumbinary 
disk was cleared: this is defined to be the point at which the 
disk inner edge /?,„ > 25AU and /?,„ > 10a. Operationally, I 
solve Equation Q] using a standard first-order explicit scheme 
on an .R'/ 2 -spaced grid (e.g., Pring le et al.ll 19861) . using 1000 
cells to span the range [0.09 AU, 2500AU] . 

3. RESULTS 

Figure Q] shows the evolution of the disk surface density in 
a characteristic model (the standard model with a = 1AU) for 
both the strong and weak wind cases. In both cases the tidal 
torque from the binary slows disk evolution, increasing the 
lifetime by a factor of ~ 2.5-3 relative to an otherwise iden- 
tical single-star disk. The most obvious difference between 
the two models is that the disk lifetime fdisk is almost an or- 
der of magnitude shorter in the strong wind case (1.68Myr 
versus 13.74Myr; the reference single-star disk lifetimes are 
0.57Myr and 5.51Myr, respectively). Note, however, that the 
absolute values of the disk lifetime are not significant, as fdisk 
depends strongly on the initial conditions (primarily the initial 
disk mass). It is more instructive to consider the evolution of 
the models as a function of the normalised time t / fdisk, as this 
highlights the relative importance of the competing physical 
processes and eliminates most of the artefacts introduced by 
the choice of initial conditions. 

If we compare these two models in this manner we see sig- 
nificant differences in their evolution. In both cases the binary 
separation (1 AU) is sufficiently small that the disk is optically 
thick interior to R c , and the initial accretion rate on to the 
binary is 3 x lO^M^yr" 1 . This is much larger than the pho- 
toevaporation rate (henceforth M w ) in the weak wind case, but 
comparable to the wind rate in the strong wind case. Conse- 
quently the two models follow qualitatively different evolu- 
tionary sequences. In the strong wind case the photoevapora- 
tive wind almost immediately overwhelms the accretion flow. 
The inner edge of the disk remains at ~ 4.2AU for much of the 
disk lifetime, but is progressively eroded by the wind. Over 
the disk lifetime most (~ 70%) of the disk mass is removed 
by the photoevporative wind, with only ~ 30% accreting on to 



Circumbinary disk dispersal 

Strong wind model, a = 1AU 



Strong photoevaporative wind 



3 




10 

R / AU 
Weak wind model, a = 1AU 



100 




R / AU 

FIG. 1 . — Surface density evolution of the standard models with a = 1AU. 
The upper panel shows the strong wind case (typical of X-ray photoevapora- 
tion), with S(i?) plotted at 0.1, 0.5, 0.7, 0.85, 0.9, 0.95 & l.Ofoi*. (Here the 
disk lifetime fcjisk = 1 -68Myr). In this case the majority of the disk is removed 
by the photoevaporative wind, and the disk is cleared rapidly once the surface 
density drops below ~ lgcm~ 2 . The lower panel shows the weak wind case 
(typical of EUV photoevaporation), with S(R) plotted at the same fractions of 
the disk lifetime (fji s k = 13.74Myr). Here accretion dominates the evolution, 
and the wind only clears the disk once it has undergone substantial viscous 
evolution. 

the binary. The high mass loss rate also has a strong impact on 
the radial profile of the disk: the surface density never reaches 
the power-law profile characteristic of accretion- or decretion- 
dominated disks. Eventually the wind prevails completely, 
and the disk is then rapidly cleared from the inside out. This 
final clearing retains the "two-time-scale " behavior charac- 
teristic of photoevapo rative disk dispersal (Clark e"et al JI200 it 
lAlexander et al.l2006bl) . but with the important caveat that ac- 
cretion and angular momentum transport play only a minor 
role in the disk's evolution. 

By contrast, in the weak wind case the initial accretion rate 
exceeds the wind rate by a factor > 100, and the disk's evolu- 
tion is governed primarily by the viscosity. At R > 10AU the 
disk is close to the £ oc R~ 3 / 2 power-law of a decretion disk 
and, as in the single-star case, only once the disk has under- 
gone substantial viscous evolution does the photoevaporative 
wind trigger disk clearing. The inner edge of the disk remains 
at a fixed position throughout, and although the final clearing 
again shows the characteristic two-time-scale behaviour, here 
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FIG. 2. — Disk lifetimes (disk plotted as a function of binary separation 
a. For clarity, all the curves are normalised to the lifetime at a = 1AU. The 
single-star models, plotted for reference, have the same viscous time-scales 
as the standard model set. The upper panel shows the strong wind models; 
the lower panel the weak wind models. The sharp increase in disk lifetimes 
for a < 0.25R C is caused by the "switch" to direct photoevaporation of the 
disk once the binary has cleared a sufficiently large inner cavity. 



the majority of the disk has been accreted on to the binary, 
with only a small fraction (7%) removed by photoevapora- 
tion. 

Further insight can be gained by looking at Figure|2] which 
shows how the disk lifetime varies as a function of binary 
separation a in the different models. The most prominent fea- 
tures are a sharp increase in the disk lifetime at small sepa- 
rations (a < 0.25/? c ), and a progressive decrease in disk life- 
times at large separations (a > 5-10AU). The increase in disk 
lifetimes at small a occurs because the photoevaporative wind 
has a characteristic radius: when the binary clears a cavity in 
the disk that is larger than the critical radius R c , the wind is 
driven by direct irradiation of the inner disk edge. This in- 
creases the efficiency of the wind, and clears the disk more 
rapidly. The increase in f^k for small a is much larger when 
accretion on to the binary is suppressed, and the critical value 
of a also increases for a lower mass ratio q. This effect is also 
more pronounced in the weak wind models. This is primar- 
ily due to the more prominent role of accretion in the weak 
wind models, but is also partly because direct irradiation in- 
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creases the weak wind rate 2 by a factor of ~ 10. By contrast, 
in the strong wind case photoevaporation dominates the evo- 
lution throughout and fdisk ~ Md/M w , so we see only small 
variations in disk lifetime at small a. 

At large a (> 5-10AU) photoevaporation also dominates 
the disk evolution: accretion in these large disks is not effi- 
cient (as the viscous time-scale is long), and M w increases for 
larger inner cavities. Consequently circumbinary disk life- 
times decrease to large a. This effect is particularly notable 
in the strong wind case: photoevaporation rapidly erodes the 
disk, and the resulting lifetimes are a factor > 10 shorter than 
corresponding single-star disks. 

4. DISCUSSION 

The 1-D models presented here are obviously simpli- 
fied, and neglect several potentially important issues. Non- 
axisymmetric features such as tidal accretion streams can only 
be modeled in an ad hoc manner, and I do not include the 
orbital modulation of the radiation fields that drive photoe- 
vaporation. This simplified model also cannot account for 
eccentricity or inclination of the binary orbit, and neglects 
both the torques and shielding effects of circumstellar disks 
around each individual star. More detailed two- and three- 
dimensional calculations are required to investigate these ef- 
fects, but these simple models still provide an important initial 
insight into how young circumbinary disks evolve. 

Perhaps the most intriguing result is that photoevapora- 
tive winds leave a potentially observable signature on the 
circumbinary disk distribution. As seen in Figure |2j all the 
models show an increase in disk lifetimes (by a factor < 2.5) 
for a < 0.257? c . Observations of circumbinary disk lifetimes 
can therefore potentially measure R c , which in turn would tell 
us the temperature (and or igin) of the photoe vaporative flow. 
Current observations (e.g.. lKraus et al.l l2012) do not yet pro- 
vide a useful census of the young binary population at ^AU 
separations, but future such surveys will offer an important 
diagnostic for protoplanetary disk dispersal. I note also that 
my results are broadly consistent with observations of debris 
disks, which are much more frequ ent around close (a < 3AU) 
binaries than at wider separations dTrilling et al.ll2007l) . 

More generally, the fact that photoevaporation dominates 
circumbinary disk evolution for high wind rates suggests that 

2 This is due to the increased efficiency of radiative transfer in the EUV 
wind once the inner disk is removed (Alexander et al. 2006a). 

3 Most interesting is perhaps DF Tau: this near-equal mass (q = 0.9) bi- 
nary has a projected separation of 10.6AU, a disk mass of 4 X 10" 4 M Q . 
and L x ~ 1 X 10 30 erg s _1 Walter & KuHTIMIjohns-Krull & Basnl[T99l. 
This implies an X-ray photoevaporation rate of 7 X lO~' , M0yr~ 1 IQwen et alj 



individual binary systems c an be used t o test photoev apora- 
tive wind models directly. iRraus et all (1201 ll 1201 2l) found 
10 Class II/III binaries in Taurus-Auriga with projected sep- 
arations less than 50AU and detected disks (either circum- 
stellar or circumbinary). 7 of these 10 objects have been 
detected in X-rays , with luminosities Lx ~ 10 29 -10 30 ergs _1 
(Giide let al.l2007l), but 9 of the lOhave e stimated disk masses 
< 10" 3 M g ( kndrews& Williamsl [2005). The existence of 
so many low-mass circumbinary disks at an age of l-2Myr 
suggests that they are not subject to strong photoevaporation, 
as the survival times for such disks are very short 3,4 . If we 
divide the observed disk masses by the stellar ages we can 
set a conservative upper limit to the photoevaporation rate of 
M w < 10~ 9 MQyr _1 . This is at least an order of magnitude 
smaller than predicted by models of X-ray photoevaporation, 
but roughly consistent with EUV photoevaporation models. 
More detailed modeling and further observations of close bi- 
naries are still required, but the existence of large numbers of 
low-mass circumbinary disks seems to place a strong upper 
limit on the efficiency of disk photoevaporation. 

Finally, it is interesting to consider the consequences of 
these results for planet formation. We now know that gi- 
ant planets around close binaries are relatively comm on 
(iDoyle et al.ll20Tlt IWelsh et al.ll20H lOrosz et alJl2012allbh . 
but little is known about their formation conditions. My re- 
sults suggest that circumbinary disks around close (a < 1 AU) 
binaries are longer-lived than otherwise identical disks around 
single stars. If disk lifetimes set a limit on the time-scale for 
planet formation, this suggests that planets can form read- 
ily in disks around close binaries. However, photoevapora- 
tion causes a decline in circumbinary disk lifetimes with in- 
creasing binary separation, implying that circumbinary plan- 
ets around wide (a > 10AU) binaries should be rarer. 
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2011) and thus a disk lifetime of ~ 5 X 10 4 yr, which is uncomfortably short 
compared to the l-2Myr age of the Taurus-Auriga association. Note, how- 
ever, that some Taurus binaries may be significantly younger than the single 



stars in the association I Kraus & Hillenbrand 2009). 

4 Also of interest is HD98800B, in the ~ lOMyr-old TW Hya associa- 
tion IFurlan et al.ll2007t lAndrews et alj|20Tol) . This a ~ 1AU binary has a 
3 X IO^Mq circumbinary disk which extends from 3.5-15AU (where it is 
tidally truncated by HD98800A), and Lx = 1.4 x 10 29 ergs~' ( Kas tner et all 
2004). Again, the disk lifetime when subject to photoevaporation at the rate 
predicted by Owen et al. (201 1), ~ 6 X 10 5 yr, is much less than the system's 
age. 
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